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A FINE POINTING SYSTEM FOR THE LARGE SPACE TELESCOPE 


INTRODUCTION 


The Large Space Telescope (LST) is to be an astronomical facility developed by 
the National Aeronautics and Space Administration (NASA) under the direction of the 
George C. Marshall Space Flight Center (MSFC), Huntsville, Alabama. The LST is the 
next step for optical stellar space astronomy after the Orbital Astronomical Observatories 
(OAO) and the Apollo Telescope Mount (ATM). It will be designed as a general-purpose 
facility, capable of utilizing a wide range of scientific instruments. The LST is expected 
to contribute significantly to studies relevant to the origin and structure of the universe, 
the study of energy processes that occur in galactic nuclei, the study of early stages of 
stellar and solar systems, and observation of such highly evolved objects as supernova 
remnants and white dwarfs. The LST will be capable of viewing galaxies 100 times fainter 
than those seen by the most powerful ground-based optical telescopes. It will weigh 
between 9000 and 11 000 kg and have a length of 12 to 16 m and a diameter of 3.6 to 
4 m. The basic LST elements (Fig. 1) are the Optical Telescope Assembly (OTA), the 
Scientific Instruments (SI), and the Support Systems Module (SSM). The most important 
optical element in the OTA will be a diffraction-limited primary mirror approximately 
3 m in diameter. The potential scientific instruments include a diffraction-limited camera 
and a low and a high dispersion spectrograph. All support systems such as control 
actuators, electronics, the power unit, and the data transmission assembly are combined 
in the SSM. Solar panels will supply the electrical power to the LST. The spacecraft will 
orbit earth at an altitude of 600 to 800 km at an inclination of 28.5 deg. Further details 
are given in Reference 1 and an LST briefing to industry.* 

One of the main problems in the early LST system analysis has been the image 
motion stabilization to within 0.005 arc sec rms necessary to achieve the maximum 
benefit from the diffraction-limited telescope. The primary stabilization approaches are 
the secondary mirror control with the spacecraft stabilized to +1 arc sec and the fine 
body pointing of the entire spacecraft to 0.005 arc sec rms. Several studies have been 
concerned with the body pointing. Johns [2] has shown that performance of about ±0.5 
arc sec is achievable using single-gimbal control moment gyros only. Hartter et al. [3] 
have studied a control system with double-gimbal control moment gyros and reaction 
wheels. Under the influence of the gravity-gradient disturbance torques, a pointing 
accuracy of ±0.001 arc sec was found. Proise [4] has investigated a control system with 
reaction wheels only. Including gravity-gradient torques and sensor noise, a minimum 
pointing error of 0.0037 arc sec rms was reported, indicating the important impact of 
sensor noise. Thus, it has been shown that performance comparable to the required 0.005 
arc sec rms might be achievable. However, critical aspects remain, e.g., gyro and reaction 
wheel vibrations,' bending modes, nonlinearities, etc., all of which may degrade the 
pointing stability. 


* LST Briefing to Industry on Phase A Results, Program Development Directorate, Marshall 
Space Flight Center, Huntsville, Ala., November 1972. 
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Figure 1 . LST spacecraft configuration. 


This report is divided into six major sections. The first section contains a 
description of the LST model and the general guidelines and assumptions. In the second 
section, the control system design philosophy is discussed, and a controller is constructed 
for the complete cancellation of the disturbance torques. The third section deals with the 
asymptotic stability and the transient analysis of the control system characterized by the 
impulse response. In the fourth section, the steady-state response of the system to noise 
excitations is investigated. The steady-state standard deviation (rms) of the LST attitude 
motion is optimized with respect to the system’s characteristic frequency featuring the 
Disturbance Accommodation Standard deviation Optimal Controller (DASOC). In the 
fifth section, the sensitivity of the DASOC caused by variations in disturbance frequency 
is found from the corresponding frequency response. Finally, in the sixth section, the 
pointing stability achievable by the DASOC is compared with the results, using a 
Proportional Controller (PC), a Proportional Integral Controller (PIC), and a Proportional 
Integral Squared Controller (PI 2 C). In addition to the factor of 3 improvement of the 
PI 2 C over the PC, the DASOC offers a possible factor of 30 improvement in pointing 
stability. Furthermore, the corresponding characteristic frequency will be far below the 
first bending mode of the LST structure. The report appendix summarizes the spectral 
density analysis and the covariance analysis of linear differential equations driven by 
white noise. In this report, the covariance analysis will be used primarily. 


DESCRIPTION OF THE LST MODEL 


The LST (Fig. 2) consists primarily of the OTA, the SI, and the SSM. The LST 
spacecraft has a nearly cylindrical, beam-like shape. The principal axes frame S , with its 
origin at the center of mass, is defined as follows: The x g -axis corresponds to the 

telescope axis, the y s -axis is located along the solar wing axis, and the z s -axis completes 

the orthogonal frame. The fine guidance system will maintain the frame S close to a 
given inertial reference frame I . Thus, the transformation matrix between frame 1 and 
frame S reads as 


IS A 


1 -i> 6 

\ P 1 -4> 


S (p 1 


<t>< 1 , 0«1 , 1 


( 1 ) 


where <t> is the roll angle, 9 the pitch angle, and \j/ the yaw angle. High accuracy is 
required in the pitch and yaw motion since these motions affect the pointing stability of 
the telescope axis. The angular velocity vector w of the LST relative to the frame I is 


u> = [<t> 6 ^ 


( 2 ) 
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Figure 2. Basic LST elements and principal axes. 


The LST will be assumed to be a rigid body; this is a legitimate assumption for 
the control system analysis because the structural frequencies are considerably higher 
than those of the control system. The reference cited in the footnote on page 1, for 
instance, indicates that the bending frequencies of the LST are all greater- than 5.7 
rad/sec. On the other hand, the characteristic, frequency of the fine guidance system will 
be less than 0.1 rad/sec. Thus, the rigid body model will be a reliable approximation. The 
inertia tensor I of the LST is given by 



where l x = 14 656 kgm 2 , I y = 91 772 kgm 2 , and I z - 95 027 kgm 2 in the on-orbit 

configuration with extended light shield and solar wings. The corresponding mass totals 
9380 kg. 

Within the SSM, there will be three symmetric reaction wheels with inertia tensors 



and the angular velocities relative to frame S: 



The reaction wheels, then, are mounted parallel to the axes x s , y s , z s , respectively. 


5 



The spacecraft’s equation of motion follows now from the Euler equations, given 
in Reference 5. 


3 

I • to + oj X 1_ ' oj + ^ [a; X coj tr 1| + 2co j X I- ■ co 

_ 1=1 


+ li f£i + ^i X Ij ' coj ] = M 


( 6 ) 


and 


li " — + — * Ij ' w + w X co j tr Ij + 2w j X I j ' co 

+ Ii ' f^i + ifi x li ‘ ^i = Mj , i = 1(1)3 , (7) 

where M is the external torque acting on the LST and Mj is the internal torque on 

the reaction wheels. Assuming small motions 4>< 1 , d<\ , \]/<\ and small reaction 

wheels Ij x ^I x , > ^3z^^z > as we ^ as tlie ^ act that th e reaction wheels have one 

degree of freedom only, the following scalar equations are obtained from the vector 
equations (1) through (7): 


r x* + 

hx^l = M x > 

(8) 

l y 0 + 

hyh = M y > 

(9) 

+ 

! 3z **3 = M z > 

(10) 

l \x u \ 

= M lx > 

(11) 

l 2y v 2 

= M 2y , 

(12) 
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and 


l 3z v 3 ~ M 3z 


(13) 


Clearly, these are three decoupled sets of equations and thus only a single axis analysis is 
necessary. Because of the high accuracy pointing requirements on the telescope axis, the 
pitch motion is chosen. 

The internal torque, M 2y , can be written as 

^2y ” ^delay + ^drive _ 


where 


M delay " ' d 0 v 2 


(15) 


is some delaying torque, e.g., viscous friction, and M drive is the driving motor 
torque, M dr j ve , which is a linear function of the control variable u and compensates for 
the delaying torque M de j a y , has the form: 


M drive = <V2 - ku 


(16) 


Thus, from equations (9), (12), and (14) through (16), it follows that 

l y 6 = ku + My (17) 

In equation (17), the external torques M y are not specified; they depend on the external 
environment in space. The maximum values of the respective external torques are given in 
Reference 1 : 

Gravity-gradient 0.2200 N-m 

Magnetic 0.0500 N-m 
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Aerodynamic 


0.0001 N-m 


Solar pressure 0.0004 N-m 

Clearly, the gravity -gradient and the magnetic torques are the significant disturbance 
torques, while the aerodynamic and the solar pressure torques are negligible. 

. The gravity-gradient torque acting on the LST spacecraft in a circular orbit with 
twice orbital rate w reads as 


Mg = e v X 1 


e v 


(18) 


where e y is the unit vector in the direction of the orbital radius vector. In the inertial 
frame O, Figure 3, e y is 

o?v = (cos i sin r sin i sin r cos t)^ (19) 


where i is the orbital inclination and t an angle locating the vehicle in orbit. The 
transformation matrix between frame O and the reference frame 1 is the constant matrix 



a ll 

a i2 

a 1 3 


a 22 

&23 

a 31 

a 32 

a 33 


( 20 ) 


Applying the transformation matrices (1) and (20) to the vertical vector (19), one obtains 
in frame S with 0^1 , 6< 1 , 0^1 : 


Sv = S 'A 


IO, 


o?v 


^ 1 « 
&12 
a 13 


cos i sin r + 
cos i sin t + 
cos i sin r + 


a 2 i sin i sin r 
a 22 sin i sin r 
a 23 sin i sin r 


+ a 31 cost 
+ a 32 cost 
+ a 33 cos r 


( 21 ) 
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Figure 3. LST orbit and inertial frames. 

Using equations (3) and (21), the gravity-gradient torque (18) can be written 

0 Z - I y )(g x l + g x 2 cos2t) 

22 , .2 

Mg = (I x - I z )(gyi + gy 2 cos 2r) (22) 

(l y - I x )(g z i + g z2 cos2r) 

where g vi < 1 , g,,; < 1 , g 7i < 1 , i = 1, 2, are introduced as convenient abbreviations. 



The magnetic torque acting on the LST is 


M m = m X B 


( 23 ) 


where jn is the LST dipole moment and J3 is the earth’s magnetic field intensity. The 
dipole moment of the LST is composed of a constant moment inherent in the structure 
and a nonconstant moment generated by desaturation magnets. In each case, the dipole 
moment will usually be fixed within the LST and may be piecewise constant in time, 
depending on the desaturation procedure. The dipole moment is characterized by a unit 
vector e m and the magnitude m 0 : 


m 


m. 


®m 


(24) 


The magnetic field intensity is, according to McElvain [6], in frame O with the 
y 0 -axis approximately corresponding to the earth’s magnetic axis: 


oB 


3B 0 


sin i cos i sin 2 t 
sin 2 i sin 2 r - 1/3 
sin i sin r cos t 


(25) 


Using the equations (1) and (20), one obtains from equations (23) through (25) the 
magnetic torque in frame S: 



h x 1 

+ 

h x2 

cos 2r 

— m = 3m o B o 

h yl 

+ 

h y2 

cos 2 t 


^zl 

+ 

^z2 

cos 2r 


where h x ^ , h y j , h z j , i = 1,2, depend on the unit vector e m , 
matrices, , and the orbital inclination i . 


(26) 


the transformation 


10 



The small aerodynamic and solar pressure torque may be modeled by a random 
torque acting on the LST: 


Mas t^asx^ ^ aS y(t) M asz (t) ^ 


(27)' 


This random torque will be a stationary white noise process with zero mean and constant 
spectral density. 

The external torque My acting on the pitch axis can now be found from 
equations (22), (26), and (27) as 


where 


and 


= ly[0 + 7 

cos (co t + X ) + 

S(t)] 

(28) 

3ce 2 *X - h 
'8 Iy 

3m 0 B 0 
8 yl + i 

y 

h yl ’ 

(29) 

3w 2 " ^z 

8 l y 

3m 0 B 0 
g y2 + i 

y 

h y2 ’ 

(30) 

+ X = 2 t 

% 


(31) 

^asyW 

V 

' (0 , q) . 


(32) 


In equation (31), the twice orbital rate oj appears again together with some 
phase angle x • The white noise process s(t) used in equation (32) is characterized by 
zero mean and spectral density q . The coefficients j3 , y , and the phase 
angle x need not be specified, because they will be generated on-line in the controller. 
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The only restriction on j8 , y , and x is that they are constant or, at most, varying from 
time to time. The twice orbital rate is approximately u> = 2 X 10' 3 rad/sec for an 
altitude of 800 km and the spectral density, q, is estimated to q = 1 X 10 " 1 2 arc 
sec 2 /sec 3 . Finally, one obtains from equations (17) and (28) the system equation, 


9 — bu + 0 + y cos (cot + x) + S 


(33) 


where 


b = k/I y 


The fine guidance sensor, part of the SI and shown in Figure 4, consists of two 
image dissectors which develop error signals as a function of the vehicle attitude motion. 
However, the sensor signals are disturbed by electronic noise generated primarily in the 
photocathode tube. The sensor signal y corresponding to the pitch motion 6 can be 
represented by 


y = 0 + v , v ~ (0 , r) 


(34) 


where the sensor noise v will be a stationary white noise process with zero mean and 
spectral density r . This spectral density r has been estimated by Kolsman Instruments, 
as shown in Reference 4. It is a function of the brightness of the guide star (m y ) and 

the amount of the guide star offset from the telescope pointing axis ( p ) . The spectral 
density can be expressed mathematically as 


r - 9.959 X lO' 18 (2.5 1)”^ p 5 ' 48 (35) 

where r is given in arc sec 2 seconds and p in arc minutes. 

For a typical guide star of magnitude m y = 12 and a guide star offset p = 20 
arc min, the spectral density turns out to be r = 8.394 X 1CT 6 arc sec 2 sec. 

Finally, the description of the LST model will be represented by the linear, 
time-invariant system, 

x = Ax + Bu + Fw + Gs (36) 
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Figure 4. Fine guidance sensor. 



and 


y = Cx + v 


(37) 


using the state vector 


x = [fl (?] T 


(38) 


and the deterministic, scalar disturbance 


w = |3 + 7 cos (cot + x) 


(39) 


The matrices A, B, C, F, G are easily obtained from equations (33), (34), (38) 
and (39): ■ 


A 


F 


0 

1 

B 

= 

_0 

0 



" 0 " 

> 

G - 

0 

_1_ 



1 


0 

b 


C = [1 0 ] , 


(40) 


A block diagram of the LST model is given in Figure 5. 


DISTURBANCE ACCOMMODATION CONTROL SYSTEM DESIGN 


For the LST modeled by the system (36) and (37), a controller must be designed. 
The design philosophy follows immediately from some principal phenomena found by 
Proise [4] . Using a lead network and a PC for the system (36) and (37), the LST 
pointing stability depends on the characteristic frequency a of the closed-loop control 
system (Fig. 6 ). With an increasing characteristic frequency, the response to the 
disturbance torques decreases whereas the response to the sensor noise increases. Thus, a 
unique optimal pointing stability can be found. For a further improvement of the 
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Figure 5. Block diagram of the LST model. 


pointing stability, either the response to the disturbance torques or the response to the 
sensor noise has to be reduced by more sophisticated controllers. Because the sensor 
noise is a stochastic process, there is very little prospect for getting a large improvement 
by a different controller. On the other hand, the disturbance torques are essentially 
known, deterministic functions which might be cancelled out by a well-matched 
controller. A reduction of the disturbance torque response offers a decrease in the 
characteristic frequency, which produces the additional advantage of reducing the 
interaction of the pointing and structural bending motion. In Reference 4, Proise has 
substituted for the PC the PIC and PI 2 C, resulting in a factor of 3 improvement in 
pointing stability and a factor of 9 reduction in optimal characteristic frequency. 
Principally, the PIC or PI 1 C, are well matched for the cancellation of constant or ramp 
disturbances, respectively, but they fit only approximately for sinusoidal disturbances. 
Therefore, a more sophisticated controller will now be designed using the accommodation 
of external disturbances in linear control systems recently developed by Johnson [7, 8, 9] 
and Davison. [10] The fundamental idea of the disturbance accommodation is the 
modeling of the disturbance function by differential equations. For this purpose, the 
shape of the waveform of the expected disturbance function must be known. This means, 
in general, that only the frequencies in the Fourier expansion of the disturbance 
function, and not the amplitudes and phase angles, give the essential and necessary 
information for the disturbance modeling. The disturbance differential equations are then 
used to complete the system differential equations, and a common observer is used for 
the estimation of state and disturbance. The state estimation is fed back by the regulator 
to obtain stability and satisfactory transient response, and the disturbance estimation is 
fed back by the accommodator to cancel out the disturbances. 
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and the disturbance model as 


i = Dz (43) 

w = Hz , (44) 


Johnson has used the following control law, 


u = u c + u r 


(45) 


with the counteraction mode, 


u c = -rw , (46) 

and the regulator mode, 

u r = -L'x . (47) 


A cancellation or counteracting control u c exists if 


F = Br 


(48) 


for some matrix V . The estimations 'X s , 'z' can be obtained from the 
observer-accommodator 


"/X 

X 


l 

- 1 
A-BL + K, C | 

0 


X 


Kx 

A 

z 


_ k 2 c f 

D 


x 


_ K 2_ 


(49) 


The controller parameters are the elements of the matrices Kj , K 2 , and L . Johnson 
recommends that a fast observer be implemented by K, and K 2 and an optimal 


17 



regulator by L . However, these recommendations cannot be used for the LST because a 
fast observer would amplify the sensor noise which was not considered in Johnson’s 
system (41) and (42). Davison’s approach to the disturbance accommodation is not as 
general as Johnson’s approach. Davison uses the complete state for the feedback in the 
regulator loop which is not available for the LST. Furthermore, Davison guarantees only 
that there will be no output response for an external disturbance input, leaving the 
possibility of an internal steady-state motion which might excite structural vibrations. 
Therefore, Johnson’s approach will be used here. 

First of all, the differential equations modeling the deterministic disturbance given 
in equation (39) must be added to the system equations (36) and (37). Since the 
disturbance is a linear combination of a constant and a sinusoidal term, the following 
disturbance state equations are used: 


Z| = 0 


(50) 


z 2 ~ 0 JZ 3 (51) 

z 3 = -wz 2 . (52) 


With these equations, the matrices in the vector equations (43) and (44) read as 


D 


' 0 0 0 ‘ 

0 0 W 

.0 -co 0 


H = [110] 


(53) 


Further, the controller parameters will be specified as: 
regulator parameter matrix, 

L = [c d] ; (54) 

observer parameter matrix, 

= [k, k 2 ] T ; (55) 

accommodator parameter matrix, 

K 2 — [k 3 k 4 k 5 ] ^ . (56) 
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Now all the matrices of the controller represented by equations (45) through (49) are 
specified. Combined with this controller, the open-loop system (36) and (37) results in 
the following closed-loop system: 


Ax 

- BLx - 

FHz 

+ Fw + Gs , 

(57) 

(A 

- BL + 

Kj C) 

x - K, Cx - K, v , 

(58) 

Dz 

+ K 2 Cx 

- k 2 

Cx - K 2 v . 

(59) 


A block diagram of the closed-loop fine pointing system is shown in Figure 7, 
representing the three different loops of the controller: 

1. The observer for reconstruction of the complete state 

2. The regulator for the state feedback 

3. The accommodator for cancellation of the disturbance torques. 

Equations (57) through (59) can be condensed by using a generalized state 
vector X ~ [x^x^z^]^, a generalized deterministic disturbance W = w, and a 
generalized noise excitation V = [sv]^ : 

X = AX + FW + GV . (60) 


The matrices A , F , and G found from equations (40) and (53) through (59) are listed 
in Table 1 . 

The initial state vector of the differential equation (60) is given by 
X(O) = X 0 (61) 

and will be introduced as a Gaussian random vector, 

X 0 ~ (X 0 , Po) , (62) 
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Figure 7. Block diagram of the fine pointing system. 
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CONTROLLER 










TABLE I . FINE POINTING SYSTEM MATRICES 
DEPENDING ON THE FREE PARAMETERS 



where X 0 is the mean vector, and P 0 is the covariance matrix. If no initial information 
on X 0 is available, X 0 = 0 and P 0 = E are useful assumptions. Here E is the unit 
matrix. The deterministic disturbance remains unchanged and follows immediately from 
equation (39) as 


W = 0 + 7 cos (cot + x) 


(63) 


The noise excitation vector is characterized by the white noise process 


V ~ (O , S) , 



(64) 


which is composed of the scalar processes (32) and (34). 

The existence of the controller represented by equations (45) through (49) 
depends on the controllability of 4he system (41), as well as on the observability of the 
system (41) through (44) and the cancellation condition (48). All three conditions will 
now be checked, using the matrices (40) and (53). The controllability matrix has the full 
rank, 


rank [B I AB] 


0 


rank 


b 


b 

0 


2 


(65) 
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proving controllability. The observability matrix corresponding to the system and 

T T T 

disturbance state vector [x 1 z ] 1 has the full rank, 


rank 



r _ 


_ 


r i 


C T 
0 _ 

• 

A FH 
0 D 

T 

C T 

_ 0 _ 


A FH 

y(n-l ) 

'c T * 

" 

O D 


0 

- 

- 


= rank 


1 

0 

0 

0 

0 


0 0 0 

1 0 0 

0 1 0 

0 1 0 

0 0 -cj 


0 

0 

0 

■to 

0 


5 


( 66 ) 


proving the observability of the system state and disturbance state. The exponent, n , is 
the order of the system defined by equations (41) and (43). The matrix r follows as 


r 


b 


(67) 


proving the counteraction mode. 

So far the controller has been designed and its existence has been proved. Seven 
parameters have been found in the matrices (54) through (56) which are free with respect 
to the disturbance accommodation mode. These parameters will be constrained by the 
stability conditions and can be optimized with respect to the transient response and/or 
the noise response. The steady-state response to the deterministic disturbance (63) is 
identically zero as a result of the disturbance accommodation system design. 


STABILITY AND TRANSIENT RESPONSE 


The asymptotic stability and the transient response of the undisturbed closed-loop 
system depend only on the homogeneous part of the differential equation (60), which is 
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( 68 ) 


X = AX , X(0) = X 0 

The asymptotic stability will be checked by applying the Lienard-Chipart stability 
conditions on the characteristic equation, 

det(XE - A) = X 7 + 3i X 6 + a 2 X s + a 3 X 4 + a 4 X 3 + a 5 X 2 + a 6 X + a 7 = 0 . (69) 
These conditions, according to Porter [11] are simply 

a 7 > 0 , a 5 > 0 , a 3 > 0 , a 3 > 0 , H 6 > 0 , H 4 > 0 , and H 2 > 0 (70) 

where Hj is the ith principal minor of the Hurwitz matrix 


ai 

1 

0 

0 

0 

0 

0 


a 3 

32 

a x 

1 

0 

0 

0 


a 5 

a 4 

a 3 

a 2 

a, 

1 

0 


^7 


3s 

a 4 

a 3 

a 2 

*1 • 

(71) 

0 

0 

37 

a 6 

35 

a 4 

a 3 


0 

0 

0 

0 

37 

a 6 

a 5 


0 

0 

0 

0 

0 

0 

a 7 



The evaluation of the characteristic equation (69) using A from Table 1 results in 

det(XE - A) = (X 2 + bdX + bc)[X 5 - k,X 4 + (w 2 - k 2 )X 3 

- (kiOj 2 + k 3 + k 4 )X 2 - u)(k 2 a> + k s )X - co 2 k 3 ] =0 . 

(72) 
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Clearly, the two eigenvalues of the system associated with the regulator and the five 
associated with the observer with accommodator are completely separated. Therefore, the 
stability conditions (70) can be applied to the second-order and the fifth-order 
characteristic equation within the seventh-order characteristic equation (72) separately. 
The results are listed in Table 2. 


TABLE 2. STABILITY CONDITIONS DEPENDING 
ON THE FREE PARAMETERS 


NO. 

STABILITY CONDITION 

1 

be > 0 

2 

bd > 0 

3 

. 

k.j < 0 

4 

k 3 <0 

5 

k 1 W 2 + k 3 + k 4 < 0 

6 

k 1 k 2 + k 3 + k 4 > 0 

7 

| (k 1 k 2 + k 3 + k 4 ) k 2 + (k^kg + k 3 ) co 2 | k 4 co 
+ | (k^ + k 3 + k 4 ) (k 3 + k 4 ) - (k 1 k 2 + k 3 - k 4 ) k.,M 2 [ k 5 
- k 1 2 kg 2 w > 0 


The seven free parameters c , d , k, , k 2 , . . . k s are now constrained because of 
the stability conditions of Table 2. But within these constraints they remain free and can 
be used to optimize the responses of the closed-loop system with respect to initial 
conditions or impulses and/or with respect to steady-state noise disturbances. Five 
different optimization approaches are listed as follows: 

1. Nonop timal transient response, wherein all seven free parameters can be used 
for optimization of the steady-state response. 

2. Suboptimal transient response with two parameters optimal in the sense of a 
quadratic performance index (Johnson’s approach of the optimal regulator). Then five 
parameters are left for the optimization of the steady-state response. 
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3. Suboptimal transient response with two additional parameters optimal in the 
sense of state estimation (optimal regulator and Kalman filter). Three parameters are left 
for the optimization of the steady-state response. 

4. Suboptimal transient response with six parameters optimal in the sense of pole 
assignment (eigenvalues proportional to a characteristic frequency). One parameter can be 
used for the optimization of the steady-state response. 

5. Optimal transient response using all seven parameters; nonoptimal steady-state 
response. 

All optimizations, except the optimal regulator, Kalman filter, and pole assignment, 
are pure parameter optimizations, which must be done by computer search methods. 

For the LST fine pointing system, the suboptimal transient response in the sense 
of pole assignment will be used so that interactions between pointing motion frequencies 
and bending frequencies will be minimized. 

Clearly, the frequencies or the negative eigenvalues must be as small as possible 
from the interaction point of view. On the other hand, the negative eigenvalues have to 
be as large as possible to furnish a good degree of stability and a short settling time. 
Thus, an optimum will be a sevenfold negative real eigenvalue. The magnitude of this 
multiple eigenvalue has to be variable for the steady-state response optimization, and this 
variable magnitude will serve as the system’s characteristic frequency. 

The control system designed as described in the section. Disturbance 
Accommodation Control System Design, is completely controllable and completely 
observable. Therefore, the eigenvalues can be chosen arbitrarily, and they will be chosen, 
optimally as a sevenfold eigenvalue 


\ = -a , i = 1(1)7 


(73) 


where a is the sytem’s characteristic frequency. The corresponding characteristic 
equation reads as 

det(XE - A) = (X + a) 7 


= (X 2 + 2aX + a 2 )(X 5 + 5aX 4 + 10a 2 X 3 + 10a 3 X 2 


+ 5a 4 X + a 5 -) = 0 


(74) 
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By comparison of the coefficients in equations (72) and (74), the seven free parameters 


are defined 

as follows: 


c = 

« 2 /b , 

(75) 

d = 

: 2a/b , 

(76) 

ki 

11 

G/l 

P 

(77) 

k 2 

- co 2 - 10a 2 , 

(78) 

k 3 

- -a 5 /o) 2 , 

(79) 

k 4 

= 5u> 2 a - 10a 3 + a s /a> 2 , 

(80) 


and 


k s = -cj 3 + lOax* 2 - 5a 4 /cj . (81) 

Using the parameters (75) through (81) and Table 1, one obtains the matrices A, F, 
G of the closed-loop system, depending on the characteristic frequency a , rewritten in 
Table 3. The asymptotic stability of the closed-loop system is now guaranteed for all 
•positive characteristic frequencies, a > 0 , following from the stability conditions in 
Table 2 or immediately from equation (73). A stability analysis is not required later. 

OPTIMIZATION OF THE NOISE RESPONSE 


The response of the closed-loop system to noise inputs depends on the stochastic 
part of differential equation (60), written as 

X = AX + GV (82) 
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with the Gaussian initial condition, 


X c ~ (X 0 , P 0 ) , 


and the white noise excitation process, 


V ~ (O , S) , S = 


q 0 


0 r 


The steady-state response to the white noise excitation process is a stochastic solution 
process which is characterized by the zero mean vector, 


E{X(t» = X = 0 , 


and the constant covariance matrix 


E |x(t)X T (t)} = P = const 
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The covariance matrix P indicates the general information in the solution process. 
For the optimization, a criterion for the performance of the LST fine pointing system 
must be introduced. For this purpose, a scalar functional of the elements of the 
covariance matrix can be used. In particular, for the LST the standard deviation of the 
pitch motion 


o e - ct xI - \/TV 


(87) 


is an appropriate criterion directly characterizing the pointing stability. The covariance 
matrix P can be found by the spectral density analysis or by the covariance analysis of 
the stochastic differential equation (82). Both methods are outlined in the appendix. The 
covariance analysis is used here because it is better suited to systems of higher dimension. 

The steady-state covariance matrix P can be computed from the algebraic matrix 
equation, 


AP + PA T + Q = 0 


( 88 ) 


where 


Q = GSG T 


(89) 


Equation (88) has a unique solution since the system (82) is asymptotically stable. 
Although the matrix equation (88) is linear in P , its solution cannot be directly found 
by simple matrix inversion. However, the expansion in linear equations is possible 
according to Chen and Shieh [12], and an explicit solution has been given by 
Muller [13]. Here, the method of Chen and Shieh is used for the numerical computation 
of P . Further, one element of P is obtained algebraically by the method of Muller. The 
symmetric matrices P and Q are written in Table 4. 

The numerical computations are all done for the twice orbital 
rate w = 2 X 10~ 3 rad/sec, representing an altitude of 800 km, and for the sensor noise 
spectral density r = 8.394 X 10' 6 arc sec 2 sec corresponding to a magnitude 12 guide 
star with a 20 arc min offset. The spectral density of the disturbance noise may be 
either q = 0 or q = 1 X 1 0" 1 2 arc sec 2 /sec 3 . The characteristic frequency is variable 
within the region IX 10~ 3 < a < 1 X 10 2 rad/sec. Some results obtained by the 
method of Chen and Shieh are given in Tables 5 and 6. The following phenomena can be 
interpreted: The solution process is almost completely correlated. Only Pi 2 may be 
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TABLE 4. COVARIANCE MATRIX P AND 
NOISE EXCITATION MATRIX Q 



r* 






-1 











p n 

P 12 

P 13 

P 14 

P 1S 

P 16 

P 17 


0 

0 

0 

0 

0 

0 

0 



sa 


P 23 

P 24 

P 25 




0 

q 

0 

0 

0 

0 

0 













EB 


P 33 

P 34 

P 35 

P 36 

P 37 


0 

Q 

k, 2 r 

k,k 2 f 

k 1 k 3 r 

k 1 k 4 r 

k,k 5 r 


p = 

EB 

CM 

P 34 

P 44 

P 45 

P 46 

p 47 

Q = 

0 

0 

k^2r 

"2*' 

k 2 k 3 r 

k 2 k 4 r 

k 2 k 5 r 



E9 




P 55 

P 56 

P 57 


0 

0 

• k n k 3 r 

k 2 k^r 

k 3 Zr 

k 3 k 4 r 

k 3 k 5 r 



E9 


P 36 

P 46 

P 56 

P 66 

P 67 


C 

0 

k,k 4 r 

k 2 k 4 r 

k 3 k 4 r 

k 4 2 ' 

k 4 k 5 r 



a 


P 37 

P 47 

P S7 

P 67 

P 77 


0 

0 

k 1 k 5 r 

k 2 k s r 

k 3 k 5 r 

k 4 k 5 r 

k 5 2 -_ 









k l = 

-5a, 

k 2 = w 2 - 10a 2 , 

k 3 = a 

5 /a, 2 1 








k 4 = 5 co 2 a - 

10a 3 + 

a 5 /w 2 

- k 5 = 

■ -a; 3 t 

lOooa 2 

- 5a 4 /td 


generally vanishing, which means that the process of the pitch angle 0 and the pitch 
velocity 0 are uncorrelated. The variance P, i or the standard deviation oq is always 

A 

less than the variance Pj 3 or the standard deviation ofi of the pitch estimation 0 . This 

means that in addition to the dynamics of the observer or filter, the dynamics of the 
LST system itself are reducing the effect of the sensor noise to the pointing stability. 
Therefore, an optimal filter does not necessarily result in a minimal steady-state response 
of the pointing stability. Further, the following standard deviations Oq can be listed as 

functions of a and q (Table 7) . The standard deviations oq for q = 0 are always less 

than in the case of q = I O' 12 arc sec 2 /sec 3 and appear to be monotonically decreasing with 
the characteristic frequency a , while for the case q = 10~ 12 arc sec 2 /sec 3 a minimum 
obviously exists. 

To obtain more information on the variance P, , or the standard deviation oq , 
respectively, this element of P will be evaluated algebraically. 

The explicit solution for the covariance matrix P is, according to Muller [13] : 


Z H k+l,l 
k=0 


2k 


2 (- 1 ) m 

m=0 


A m Q A 2k 


-m 


P 


2 det H 


(90) 
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TABLE 5. COVARIANCE MATRICES FOR VARIOUS CHARACTERISTIC 

FREQUENCIES (q = 0) 


a = 2-10'^ rad/sec 


SOLUTION MATRIX P 


q - 0 

.256755 007 .000000 

.370910012 

.211162-007 

—.923340-010 

.756772-007 

.477409 010 
.114892-013 
.125910-009 
.530596-012 

.115417-013 

—.356745-016 

.451177-013 

.100728-015 

.293790-019 

.839400-013 

-.545610-016 

.186766-012 

.507837-015 

.125910019 

.671520-018 

—.<*40685-013 

.310578-015 

-.209850-012 

-.344154-015 

-.125910-018 

-.503640-018 

.671520-018 

a - 2-10 rad/sec 


SOLUTION MATRIX P 


q = 0 

.638533 006 .000000 

.787404-009 

.581504-006 

-.137679-007 

.115585-005 

.796502-008 
.319319-009 
.154641-007 
.771641 009 

.281528-002 

-.464230-010 

.617288-008 

.123713-009 

.366191-010 

-.257908-008 

.391574-010 

-.555768-008 

-.113013-009 

-.330225-010 

.297919-010 

.132326-008 

-.279059-010 

.305108-008 

.577609010 

.179693-010 

-.161814-010 

.886027-011 

a = 2-10’^ rad/sec 


SOLUTION MATRIX P 


q = 0 

.645474-005 .423516-021 

.792947-006 

.587276-005 

-.138201-005 

.116190-004 

.800018-006 

.322124-006 

.154525-005 

.774512-006 

.283929-003 
—.465591-004 
.619689-003 
.123940 003 
.367227-001 

—.283691-003 

.464855-004 

-.619069-003 

-.123833-003 

-.366865-001 

.366503-001 

.135063-004 
-.289487-005 
.312081 004 
,■587480-005 
.183579-002 
-.183396-002 
.922962-004 




TABLE 6. COVARIANCE MATRICES FOR VARIOUS CHARACTERISTIC 
FREQUENCIES (q = 1 X 10' 12 arc sec 2 /sec 3 ) 


Ct = 2-10 ^ rad/sec 


SOLUTION MATRIX P 

Q = 1-10-12 

2 3 

arc sec /sec 

.295924-003 .000000 

.465215 009 

.282248-003 

-.235299-007 

.270583-003 

-.113234 006 
.183709609 
-.117062-006 
.152874-009 

.101574609 
-.109411 012 
.101608-009 
-.936493-013 
.218779-015 

-.312416609 

-.313046613 

-.296688609 

.156758612 

-.623741-016 

.100067-014 

-.203169609 

-.187189612 

-.187710609 

-.628442613 

.623741616 

.249496615 

.100067-014 

Ct 2-10‘^ rad/sec 


SOLUTION MATRIX P 

q = 1-10' 12 

arc sec /sec 

.774509006 .000000 

,870835-009 

.712343-006 

—.141899607 

.128496-005 

.797332-008 
.355466-009 
.152917-007 
.791422 009 

.331474608 

-.497326610 

.667234-008 

.122333-009 

.388875-010 

-.304778008 

.416903-010 

-.802355-008 

-.111947-009 

-.350979610 

.317018610 

.154694-008 
-.304650 010 
.327994-008 
.566720610 
.190986610 
-.171983610 
.945471611 

Ct = 2-1 O' 1 rad/sec 


SOLUTION MATRIX P 

q - MO' 12 

arc sec^/sec^ 

.645534-004 .277556-016 

.793003-003 

.587333004 

-.138206-003 

.116196-003 

.800053-004 

.322152603 

.154928603 

.774641603 

.283963+002 

-.465604+002 

.619713+002 

.123943+003 

.367237+008 

-.283951+002 

.465597+002 

-.619707+002 

-.123942+003 

-.367234+008 

.367230+008 

.135091+000 

-.289592+000 

.312151+000 

.587579+000 

.183618+006 

-.183616+006 

.923336+003 




TABLE 7. STANDARD DEVIATIONS 


a e 

(arc sec) 

q = 0 

q = 10-' 2 
arc sec 2 /sec 3 

a = 2 X 10 3 rad/$ec 

1.599 X 10" 

1.740 X 10” 2 

a = 2 X 10" 2 rad/sec 

7.997 X 10" 

8.801 X 10" 

Q = 2 X 10 _1 rad/sec 

2.540 X 10' 3 

2.541 X 10~ 3 


where H is the Hurwitz matrix (71), 1 ^ the cofactor of the element h k+Ll ofH ’ 

A m = AA m _j + a m E is a 7 X 7-matrix, and a m is the corresponding coefficient of the 
characteristic equation (69). In particular, for Pi j it yields 


Pu 2 det H l^ u Q n + H 2 i [2(A 2 Q)n “ (Ai QA t ) n 


+ H 3 i [2(A 4 Q) n - 2(A, QA 3 ) n + (A 2 QA 2 )ii] 


+ H 41 [2(A 6 Q) 1 1 - 2(Aj QA S )ji + 2(A 2 QA 4 ) u - (A 3 QA3)u] 


+ H sl [2(A i QAj) 11 - 2 (A,Qa7)„ + (A 4 QaJ) u ] 


+ H 61 [2(A 4 QAj) n - (A s Qa7)i, 1 


+ H71 (A 6 QAg )u 


(91) 


All cofactors, matrix elements, and the Hurwitz determinant are listed in Table 8 
as a function of a , w , q , and r . From Table 8, one immediately obtains with 
equation (91) the following result: 
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TABLE 8. DETERMINANT, COFACTORS AND MATRIX ELEMENTS 
FOR CALCULATION OF VARIANCE P n 


DETERMINANT 

det H = 2 097 152 a 28 

COFACTORS 

H 11 = 437 088a 27 


H 12 = -43 008a 25 


H 13 * 14 336a 23 


h 14 = -10 240a 21 


H 15 = 14 336a 19 


h 16 = -43 008a 17 


H 17 = 437 088a 15 

MATRIX ELEMENTS 

Q 11 

(A 2 Q> n = (A 4 Q» 11 = lAgQ)^ - 0 


' a ,qa;i„ 

= q 


(A-|QA 3 >ii 

= 21a 2 q 


(A^QAg)n 

= -w 4 q + 21a) 2 a 2 q 

- 

(A 2 QA 2 

= 49a 2 q + 49co 4 a 2 r - 490oj 2 a 4 r + 1225a 6 r 


(a 2 qaJ> 11 

* 49co 2 a 2 q - 147cj 2 a 6 r + 735a 8 r 


(A 2 QAg) 1 1 

■ — 7tJ 2 a 8 r + 35a 10 r 


< A 3 QA 3>11 

= 441a 4 q + oA - 42to 6 a 2 + 511« 4 a 2 r 



-1470cj 2 a 6 r + 1225a 8 r 


(a 3 qa|) 11 

= — 21co 4 a 2 q + 441oj 2 a 4 q + 7a> 4 a 6 r 



— 147a; 2 a 8 r + 245a 10 r 


(a 4 qaJ) 11 

= 49co 4 a 2 q + 441a 10 r 


< A 4 QA 6>11 

= 21a 12 r 


<A 5 QA 5*11 

= co 8 q - 42aj 6 a 2 q + 441to 4 a 4 q + 49a 12 r 


< A 6 QA 6>11 

= a 14 r 
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p n = (2275 a' 3 - 6454 oj 2 a' 5 + 9888 cj 4 a' 7 - 882cu 6 a‘ 9 

+ 21co 8 a' n )q + ^-g (7875 a - 7350 w 2 a~' + 2800 cj 4 a' 3 

- 210ca 6 a' s + 5ca 8 a' 7 )r . - (92) 

Thus, it results that the variance Pi , is a linear combination of the disturbance noise 
term characterized by q and the sensor noise term given by r . For low characteristic 
frequencies, one obtains 


Pn 


21 

2048 


w 8 *- 11 


2048 


w e cT 7 r 


for a < co 


(93) 


whereas for high characteristic frequencies it yields 


Pn 


7875 

2048 


ar 


for a > co 


(94) 


The limit cases, a -*■ 0 and a -*■ °° , result both in P -*■ 00 even if the disturbance 
noise vanishes, i.e., q = 0 . Thus, there will be at least one minimum, depending 
generally on the ratio q/r . But the most interesting case is q = 0 , for which Table 7 
indicated a monotonically decreasing function of the characteristic frequency a . The 
extrema of P, i can be found from 


dP,,/da = 0 


(95) 


Using equation (92) with q = 0 , it follows from equation (95) that 

225 + 210 oj 2 a“ 2 - 240o J 4 cf 4 + 30<u 6 cT 6 - a/cf 8 = 0 . (96) 

According to Descartes’ rule, equation (96) has one or three positive roots for the 
variable coa 1 . It turns out that three positive roots exist, which are listed in Table 9. 
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The roots Nos. 1 and 3 represent the 
minima, and No. 3 is the absolute minima. 
In Figure 8, the standard deviation ag or 

pointing stability is shown as a function of 
the characteristic frequency. The absolute 
minimum found for q = 0 occurs at 
a = 1.5 X 10‘ 3 rad/sec and for q = 10' 12 
arc sec 2 /sec 3 at a = 1.8 X 10~ 2 rad/sec. 
With these optimizations, the controller is 
completely defined. The characteristic fre- 
quency of the pointing system is, in each 
case, far below the bending frequencies. 

Finally, some comments on the disturbance accommodation controller may be 
useful. The original approach of the disturbance accommodation presented by Johnson 
[7,8] is very well suited for a completely deterministic system given by equations (41) 
and (42). However, if there is some sensor noise within the system, the advantage of the 
disturbance accommodation may completely fade away because of the system’s response 
to the sensor noise. Therefore, the disturbance accommodation controller must be 
optimized with respect to the standard deviation of the steady-state noise response. Such 
controllers will be called Disturbance Accommodation Standard deviation Optimal 
Controllers (DASOC). The recommendation of Johnson to use fast settling times for the 
observer-accommodator will result in large negative eigenvalues representing high 
frequencies in the closed-loop system. However, high frequencies amplify the sensor noise 
heavily. On the other hand, the well-known technique using low frequencies for the 
reduction of the system’s sensor noise response will not work either. Even if all 
disturbances are exactly cancelled out, the sensor noise may be strongly amplified within 
the accommodator loop of the controller. Somewhere between low and high frequencies, 
there will be an optimum. As pointed out in the section, Stability and Transient 
Response, at least one free parameter of the closed-loop system must be available for the 
variation of the system’s frequencies or eigenvalues to obtain an optimum. 


TABLE 9. EXTREMA OF 
STANDARD DEVIATION 




°9 

Root No. 

a/to 

(arc sec) 

1 

0.258199 

3.903 X 10" 4 

2 

0.268824 

3.905 X 10"* 

3 

0.750895 

1.247 X 10"* 


SENSITIVITY ANALYSIS OF THE DASOC SYSTEM 


In the previous section, the DASOC was completely defined, and the steady-state 
response to constant and sinusoidal disturbance torques was designed to be zero if the 
disturbance frequency is exactly known. However, the disturbance frequency or twice 
orbital rate is slowly time-varying and the question arises as to the sensitivity of the 
DASOC to frequency deviations. The sensitivity can be determined by means of the 
closed-loop system’s steady-state response to the deterministic disturbances. 
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The response to the constant and sinusoidal disturbances depends on the 
deterministic part of differential equation (60) written as 


X = AX + FW 


(97) 


with the scalar excitation 


W = 0 + 7 cos (ft t + x) 


(98) 


where ft is an arbitrary frequency. Since equation (97) is linear, the constant and the 
sinusoidal disturbances can be treated separately. 

The steady-state response of the asymptotically stable system (97) to a constant 
disturbance 


W0 = 0 


(99) 


is also a constant response 


Xjj - const , X.0 = 0 


( 100 ) 


The substitution of equations (99) and (100) in the system equation (97) results in 

AX^ + Fj3 = 0 , (101) 


with the solution, 


X^ = A' 1 F0 = [0 0 0 0 (3 0 0] T 


( 102 ) 


The solution ( 1 02) has been obtained using the matrices of Table 3. The steady-state 
response of the pitch motion is exactly zero, 6 ^ = 0 , independent of the twice orbital 
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rate to . The estimation of the constant disturbance in the steady state is exactly correct, 
i.e., z, = 0 . 


For the computation of the frequency response of the system (97), it is useful to 
introduce complex functions. The sinusoidal disturbances will be represented by 

= 7COs(S2t -x) = y c cos £2 1 + 7 S sin f2t 


= fe im + fe' int 


(103) 


where 


f = J (7 C - ir s ) = ^ (tcosx - by sin x) - (104) 

rhe steady-state response is also a sinusoidal vector function: 

Xy = g e int + g e' lfit (105) 

where 

x 7 i = cos(f2t + Xi ) = S ci cosOt + | si sini2t (106) 

and v 

g i = i ^ci - ^Si^ » 1 = 1(1)7 . (107) 


With equations (103) and (105), it follows from the system equation (97) that 

g = (i£2E - A)' 1 F - f = Rf . (108) 
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The complex vector 


t» _ adj (i£2E - A) • F 
det(ifiE-A) 


(109) 


has elements of the form 


Ri 


x i + »yj 

uo + ivo 


( 110 ) 


Therefore, the amplitudes of the solution vector elements (106) are 


«i = 


fx? + yj 

u 0 2 + Vo 2 


( 111 ) 


Thus, the problem is reduced to the calculation of the vector. R by means of the 
matrices in Table 3. The result is presented in Table 10. For arbitrary disturbance 
frequencies il , the response is not vanishing. However, disturbances with the design 
' frequency w - £2 are cancelled out: 


X T (J2=co) =[0 0 0 0 0 y cos (cot - x) j ~ T sin(u>t - x)] ^ • (112) 


Notice that the observer-accommodator creates a correct estimation of the sinusoidal 
disturbance % = y cos(cat - x) , not only in amplitude but also in phase angle. 

The pointing" stability of the LST is a function of the pitch motion. Therefore, 
the pitch motion amplitude, = U will be evaluated in more detail now. One obtains 

from Table 10, / 


Xi = 7ce£2 4 - 7 w 2 aS2‘ 


y, = f2 s - 21 a 2 ft 3 - (ta 4 - 21 w 2 a 2 )ft , 


(113) 

(114) 
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TABLE 10. COMPLEX FREQUENCY RESPONSE COMPONENTS 


COMPONENT NO. 


COMPLEX COMPONENT 


+ ^ [(-tf-to 2 + 21a 2 ) + i7 afil 
det {ift E - A) L J 


~^ 2 (~n 2 + co 2 ) ^_ n 2_ OJ 2 + 21a 2, + i7a ^J 

det (iftE - A) 


ifl {-^ 2 + gj 2 )|~ (20a 2 _ a ; 2j + j5an J 

det <iflE - A) 


in <-n 2 + oj 2 ) r _ 5a 2 + j(cj 2 _ 10a 2)n] 

det (i^E - A) L J 


— (-S2 2 + CJ 2 ) r { _ a 2 + a 2) + i2an 1 
det (iftE - A) L 


k d — i^coke Poo ~1 

1 ! I M7 2 + a 2 + i2aft)J 

det (i£2E — A) 


kc + ifi co kj Poo ”1 

* 1 [i-n 2 + a 2 + i2afi>J 

det (i£2E - A) 


- 5co 2 a — 10a 3 + a 5 /<o 2 , k 5 - — to 3 + lOcoa 2 - 


det (ifiE - A) = [-7a£2 6 + 35a 3 fl 4 - 21a 5 £2 2 + a 7 ] 

+ i f-ft 7 + 2la 2 I2 5 - 35a 4 fl 3 + 7a 










Uo = -lan 6 + 35 a 3 ft 4 - 21 a 5 ft 2 + a 7 , (115) 

and 

v 0 = -ft 7 + 21 a 2 ft 5 -■ 35a 4 ft 3 + 7 a 6 ft . (116) 


Introducing the equations (113) through (116) in equation (111), the amplitude 
ratio £gjy of the pitch motion frequency response can be easily computed. In Figures 9 

through 11, the ratio $g/y is shown as a function of the actual disturbance 

frequency ft for various characteristic frequencies a . The disturbance accommodation 
point ft = cc is obvious. Further, one realizes that the sensitivity to frequency 
deviations is rapidly increasing with decreasing characteristic frequencies, For small 
frequency deviations 


ft - co(l + e) , e < 1 , 


the amplitude is directly proportional to the error e , having the form 


h = 


1 / « 1 2 + * 2 3 
CO 2 V K3 2 + Ki, 2 


ye 


where 


Ki = 14acj"' , k 2 = 4 - 42a 2 cj~ 2 
k 3 = -7 aw' 1 + 35 a 3 o; 3 - 21a s cu’ 5 + a 7 aT 7 , 
= -1 + 2la 2 a>' 2 - 35a 4 co~ 4 + 7 a 6 u>~ 6 


(117) 


(118) 


(119) 


To get an idea of the real magnitude of the amplitude ^ affecting the pointing 

stability, let us assume the worst gravity-gradient torque and a vanishing magnetic torque. 
For I x < I y = I z , g y 2 = 1 and m 0 = 0 , it follows from equation (30) that 
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DISTURBANCE FREQUENCY Si (rad/sec) 

Figure 10. Frequency response for various characteristic 

frequencies (2 X 10~ 2 < a < 10' 1 ). 






( 120 ) 


y 



Introducing equation (120) in equation (118) and using the identity 1 rad = 206 264 
arc sec, a sensitivity factor is defined as 




= 77349 



( 121 ) 


where $$ is given in arc seconds. In particular, for a > co , equation (122) becomes 


f fl /e = 3 248 658 w 5 a" 5 


( 122 ). 


The sensitivity factor fy'/e is plotted in Figure 12 as a function of the characteristic 

frequency. It is obvious again that the sensitivity is strongly decreasing with the 
increasing characteristic frequency a . 

The frequency response characterized by can be added to the noise 
response o e to obtain the resulting root mean square error or pointing stability: 


rmStf =' >/ V + 


(123) 


In Figure 13, the pointing stability (123) is shown for various errors of the 
disturbance frequency and vanishing noise disturbance, q = 0 . The optimal 
performances are listed in Table 1 1 . 

Thus the twice orbital rate co must be adjusted accurately within the DASOC to 
obtain full advantage. This can be realized by updating during the mission. But even if 
the orbital rate is not adjusted, resulting in e = 10' 1 , the DASOC will still meet the 
requirement on the LST pointing stability with some margin. 


COMPARISON WITH CONVENTIONAL SYSTEMS 


Thus far the DASOC has shown good behavior, but unquestionably this 
sophisticated controller results in more complex hardware and increasing costs. By 
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Figure 13. Pointing stability for various frequency errors. 



comparison with the simple conventional 
controllers, it can be determined if the 
DASOC is worth the additional expense. 
The conventional controllers PC, PIC, and 
PI 2 C designed by Proise [4] will be used 
in the comparison. Control systems with 
different kinds of controllers are compa- 
rable if identical inputs are used and the 
same outputs are observed. This means 
that according to Reference 5, aerody- 
namic and solar pressure, torques are not 
included, q = 0 , and a correct orbital 
rate adjustment within the DASOC is 
assumed, e = 0 . The standard devia- 
tion oq is equivalent to the 

rms0 obtained in Reference 5. Both quantities describe the pointing stability as a 
function of a corresponding characteristic frequency a . 

The comparison is shown in Figure 14, together with the structural bending 
frequencies given in the reference cited in the footnote, p. 1 . As mentioned earlier, the 
PPC results in a factor of 3 improvement in pointing stability and a factor of 9 
reduction of the optimal characteristic frequency with respect to the PC. Compared with 
the PI 2 C, the DASOC offers a factor of 30 improvement in pointing stability and a 
factor of 500 reduction of optimal characteristic frequency. 

The actual requirement on the LST pointing accuracy is stabilization within 0.005 
arc sec rms. Figure 14 indicates that the PC does not meet this requirement. The PIC 
meets the 0.005 arc sec with little margin, and the PI 2 C is scarcely better. On the other 
hand, the DASOC easily meets the requirement, and a large safety factor remains for the 
compensation of possible deteriorations. Some additional errors may be caused by 
aerodynamic and solar pressure torques and incorrect orbital rate adjustments. Also, the 
extreme accuracy of the DASOC allows greater tolerance on other sources of errors such 
as the Control Moment Gyro (CMG) and reaction wheel vibrations, bending modes, 
nonlinearities, etc. 

The characteristic frequency a is generally not identical with the frequencies in 
the transient response. These depend on the eigenvalues of the closed-loop pointing 
system. Therefore, it is useful to compare the eigenvalue distribution of different 
controllers. The eigenvalues can be found from the characteristic equation. According to 
Reference 5, the characteristic equation of the closed-loop system with the PC reads as 


TABLE 11. OPTIMA OF 
POINTING STABILITY 


€ 

rm$a •„ 
o mm 

(arc sec) 

icr 10 

1.27 X Kf 4 

10" 7 

5.45 X itr 4 

icr 4 

1.05 X 10 3 

10' 1 

1.99 X 10 3 


X 3 + 5.62aX 2 + 10a 2 X + 5.62a 3 = 0 


(124) 
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Figure 14. Pointing stability for different controllers. 





The corresponding eigenvalues are 


Xi = '1.18a , X 2 = -1.81a , X 3 = -2.63 a . (125) 

With the PIC, one obtains 

X 4 + 5.62 aX 3 + 1 0 a 2 X 2 + 6.75 a 3 X + 0.635 a 4 = 0 , (126) 

and the eigenvalues 

X, = -0.111a , X 2 = -3.03 a , X 3j4 = -1.25 a ± 0.57 ai . (127) 

With the PI 2 C, it follows that 

X s + 5.62 aX 4 + 10a 2 X 3 + 6.75 a 3 X 2 + 0.663 a 4 X + 0.0178a 5 = 0 (128) 

and the eigenvalues 

X I;2 = -0.056a , X 3 = -3.03 a , X 4>5 = -1.25 a ± 0.57 ai (129) 

The characteristic equation of the system with the DASOC has been developed in the 
section, Stability and Transient Response. The characteristic equation (74) is simply 

(X + a) 7 - 0 ( 130 > 

and the sevenfold eigenvalue is 

Xj = -a , i = 1(1)7 . 031) 

All eigenvalues are presented in Figure 1 5. 
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Figure 15. Eigenvalues for different controllers. 

It is clear that especially for the system with the PI 2 C, the eigenvalues are widely 
spread. For the characteristic frequency a the eigenvalues for the PI 2 C are contained in 
the region, 0.052a < X < 3a . The eigenvalue, X = 0.05a , results in an unnecessarily 
large settling time, whereas the eigenvalue X = 3a enhances the interaction with the 
structural bending modes. From these points of view, the DASOC has the optimal 
eigenvalue distribution. 

The optimal characteristic frequency of the fine pointing system with the DASOC 
has been found as a = 1.5 X 10" 3 rad/sec in the limit case q = 0 , e = 0 . This means 
a- large settling time of the order of magnitude of the orbit period. Such a large settling 
time may be acceptable during fine pointing operations, but it is not acceptable for LST 
maneuvering. However, this should not present a problem because maneuvering will be 
accomplished by CMG’s and not the reaction wheels. 
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CONCLUSIONS 


The LST fine body pointing to within 0.005 arc sec rms is the subject of this 
report. This difficult requirement on pointing stability is necessary to take full advantage 
of all the benefits from the diffraction-limited telescope. 

The LST is modeled as a rigid body with three orthogonally mounted reaction 
wheel actuators. Nonlinear coupling phenomena do not appear, and, therefore, only a 
single axis analysis is performed. The principal disturbance torques result from the 
gravity-gradient and the magnetic desaturation system. Both are characterized by a 
constant term and a sinusoidal function of time with twice orbital rate. The fine guidance 
sensor providing the measurement of the spacecraft attitude motion generates a signal 
cluttered with electronic noise. Altogether, the LST is represented by a linear 
time-invariant system with deterministic and stochastic disturbances. 

From earlier studies, it is known that the pointing stability of the LST with 
conventional controllers is limited by the external torques as well as by the sensor noise. 
Therefore, a more sophisticated controller is designed, using the recently developed 
technique of the accommodation of external disturbances to avoid the stability limitation 
caused by the external torques. Such a controller is constructed of three different loops: 
(l)an observer for the reconstruction of the complete state, (2) a regulator for the state 
feedback, and (3) an accommodator for the cancellation of the disturbance torques. The 
existence of a disturbance accommodation controller for the LST is proven by 
controllability and observability analyses. However, seven parameters or control gains 
remain to be determined. The stability analysis yields constraints on these free 
parameters, and within these constraints the parameters are optimized with respect to the 
transient response and the steady-state noise response. The transient response is found to 
be optimal for a sevenfold, real, negative eigenvalue. The magnitude of the multiple 
eigenvalue, the characteristic frequency, is then used for the steady-state noise response 
optimization. The steady-state response to the noise inputs is characterized by the 
covariance matrix and has been computed by the covariance equation widely used in 
estimation theory. In particular, for the standard deviation of the pointing motion, an 
algebraic formula was obtained. Then the optimal standard deviation was easily found. 
The result is the DASOC. In contrast to the original disturbance accommodation 
technique, the DASOC takes care not only of the external disturbances but also of the 
internal sensor noise. 

Thus far the LST spacecraft with the DASOC has the following properties: 

1. The asymptotic stability of the system is a priori guaranteed. 

2. The system has only one characteristic frequency (multiple root) far below 
the bending frequencies of the LST structure. 
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.v The constant and the sinusoidal disturbance torques with twice orbital rate are 
completely cancelled out. 

4. The steady-state standard deviation of the attitude motion or pointing 
siahility is a function of the characteristic frequency. 

5. An optimal pointing stability is maintained. 

The frequency response to an external disturbance with variable frequency is used 
to define the sensitivity of the DASOC to the orbital rate adjustment. It turns out that 
updating of the orbital rate simulated in the DASOC is necessary to obtain good results. 

In a comparison with the best conventional controller found in earlier studies, the 
DASOC offers up to a factor of 30 improvement in pointing stability, resulting in an 
optimal performance of nearly 0.0001 arc sec rms. In addition, the optimal characteristic 
frequency is reduced . by a factor of 500. This means the DASOC easily meets the LST 
pointing stability requirement of 0.005 arc sec rms, even if aerodynamic and solar 
pressure torques as well as incorrect orbital rate adjustments are admitted. Excitation of 
vibrations in the LST structure by the DASOC is not expected to be a problem. 

Further effort is necessary to apply the DASOC to more realistic spacecraft 
models and actuator systems and to perform a three-axis simulation. But the large-scale 
improvement with the DASOC allows a significant margin before the requirements of the 
LST are exceeded. 


George C. Marshall Space Flight Center 

National Aeronautics and Space Administration 

Marshall Space Flight Center, Alabama, 35812, April 4, 1973 
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APPENDIX 


RESPONSE OF LINEAR DIFFERENTIAL EQUATIONS TO WHITE NOISE 


Stochastic differential equations are widely "used in control theory. In the 1940’s, 
the analysis of linear differential equations driven by white noise was conducted primarily 
in the frequency domain, using the spectral density. Since the 1960’s, stochastic 
differential equations have been investigated in the time domain especially in connection 
with the development of the Kalman-Bucy filter using covariance matrices. Both methods, 
the spectral density analysis and the covariance analysis, will be reviewed briefly with 
respect to the application for the LST, Reference is made to James, Nichols, and Phillips 
[14] and to Bucy and Joseph [15]. 

The linear differential equation system 


x = Ax + Bw , x(t 0 ) = x 0 ~ (m 0 ,Po) 


(A-l) 


may be given by the n X 1 -state vector x(t) , the r X 1 -excitation vector w(t), 
the n X n-matrix A, and the n X r-matrix B. Both matrices will be time invariant, and 
the matrix A will have eigenvalues with negative real parts only; i.e., the system (A-I) will 
be asymptotically stable. The initial condition x(t 0 ) is given by a Gaussian random 
vector x 0 with mean value m 0 and covariance matrix P 0 - 

The excitation vector w(t) represents a stationary white noise process. The white 
noise process is Gaussian, it is a Markov process, and its integration yields a Wiener 
process. The mean vector of the white noise process vanishes, 

E jw(t)J = 0 , ' (A-2) 


and the covariance matrix function is given by 

E |w(t)w^ (r ) | = Q6 (t - t) , Q = Q T > 0 (A-3) 

The spectral density matrix of the stationary white noise process is constant, 


S w (oj) = Q = const 


(A-4) 


preceding page 


PLANK NOT FILMED ' 
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and the covariance matrix is infinite 


E |w(t)w^ (t) | = Q5(t ,t) -*■ 00 


(A-5) 


This means that the white noise process does not exist in the physical world, but it is an 
admissible idealization for system analysis similar to the application of impulses in the 
analysis of deterministic differential equations. 

The state vector x(t) characterizes a stochastic process, too, because of the white 
noise excitation of the system. Tire solution process is a nonstationary, Gaussian vector 
process and a Markov process. This means that the solution process can be described 
by the mean vector, 


E | x(t) | = $>(t , t 0 )m 0 


(A-6) 


and the covariance matrix function, 


E |x(t) x T (t)j = 4>(t,t 0 )P 0 <h T (t,to) 

t 

+ f <b(t , r) BQB T <t> T (t , r) dr , (A-7) 

to 

where <t>(t, t 0 ) = is the fundamental matrix of A. The mean vector 

E{x(t)} = m(t) and the covariance matrix Ejx(t) x^(t) J = P(t) satisfy the differential 
equations, 

m - Am , m(t 0 ) = m 0 , (A-8) 

P = AP + PA T + BQB T , and P(t 0 ) = P 0 . (A-9) 


Since the system (A-l) is assumed asymptotically stable, the steady-state solutions of the 
differential equations (A-8) and (A-9) are 
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OQ 


(A-10) 


m 


= / 

0 


e As BQB t e 


A T s ds 


Thus, the steady-state response of the system (A-l) is characterized by the steady-state 
covariance matrix , which can be found by the spectral density analysis or the 
covariance analysis. 

Used in the spectral density analysis is the fact that the steady-state covariance 
matrix can be found as an infinite integral over the spectral density matrix S x (co) of 

the solution process, 



S x (to ) dco 


(A-l 1) 


The spectral density matrix S x (tu) can be determined simply by matrix multiplications 
S x (w) = F(-oj) Q F t («) , (A-l 2) 


where 


F(cu) = (koE - A)' 1 B 


(A-l 3) 


Thus, the spectral density analysis is convenient up to the solution of the integral (A-ll). 
However, integral tables .are available for the evaluation of the diagonal elements P- of 

the covariance matrix, sometimes called variances or squared standard deviations. The 
variances Pjj can always be presented in the following form: 


‘JJ 


J_ r «„(<■>> 

2ni h n (cu) h n (-o)) 


1(1)7 


(A-l 4) 
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where g n (w) and h n (w) are the polynomials 


g n (cj) = b G w 2n " 2 + bj w 2n_4 + ... + b n-1 


(A-15) 


and 


h n (w) = a 0 co n + aj oj n 1 + . . . + a n = det(iwE - A) (A-16) 

A table for the integrals of the form (A-14) is given by James, Nichols, and Phillips [14] 
for systems of order n < 7 . 

The covariance analysis uses the linear algebraic covariance matrix equation, 

AP m + APj + BQB T = 0 , (A-17) 


following from the differential equation (A-9) in the steady state. The following theorem 
on the uniqueness of the solution of equation (A-17) is given by Lancaster [16] . 

Tire matrix equation (A-17) has a unique solution if and only if the n 2 numbers 
Xj + Vj * 0 , i - l(l)n , j = 1(1 )n , (A-18) 

where X x , X 2 , . . . X n are the eigenvalues of the matrix A. Since it has been assumed 
that ReXj < 0 for all i = 1(1 )n , it follows necessarily that Xj + Xj < 0, and, 
therefore, equation (A-17) has a unique solution. 

Although the matrix equation (A-17) is linear in , its solution cannot be found 
directly by matrix inversion. However, the expansion in linear equations is possible, 
according to Chen and Shieh [ 12] : 

2Cp + q = 0 . (A-l 9) 
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The vector p contains the n(n + l)/2 essential elements of the symmetric 
matrix , the vector q includes the n(n + 1 )/ 2 essential elements of the symmetric 

matrix BQB^ , and C is the n(n + l)/2 X n(n + l)/2 Chen-Shieh matrix. Furthermore, 
an explicit solution of equation (A-17) has been found by Muller [ 13] : 


, n-1 2k 

p - = VW Z n H k + l,l E „ <-D m A m BQ B T aT , (A-20) 

u k=0 m=0 

where A m = AA m _j + a m E is an n X n-matrix, a m = -tr(AA m _j)/m is the mth 

coefficient of the characteristic equation, H is the n X n-Hurwitz matrix, 
a nd Hj, + j j is the cofactor of the element h^j j of the Hurwitz matrix. 

Both methods, the spectral density analysis and the covariance analysis, can be 
used to compute the steady-state noise response. However, for systems of higher 
order, n > 5 , the covariance analysis makes the computer Programing easier because 
straightforward matrix techniques are available. Therefore, the LST fine pointing system 
of order n = 7 has been investigated by the covariance analysis. 
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